Colliding particles in highly turbulent flows 
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Abstract 



We discuss relative velocities and the collision rate of small particles suspended in a highly 
turbulent fluid. In the limit where the viscous damping is very weak, we estimate the relative 
velocities using the Kolmogorov cascade principle. 

PACS numbers: 

05.20.Dd Kinetic theory 
45.50.Tn Collisions 
47.27.-i Turbulent flows 
47.57.E- Suspensions 

Introduction. This paper considers collisions of small particles suspended in a highly turbu- 
lent gas. The collisions of these particles can facilitate aggregation of the suspended particles. 
This process may be relevant to the precipitation of rain from turbulent cumulus clouds p.] , and 
to the formation of planets by aggregation of dust particles suspended in the gas surrounding a 
growing star [2]. 



The suspended particles are characterised by a dimensionless measure of the importance of 
inertia, termed the Stokes number: St = I/7T, where r is a correlation time of the flow and 7 is 
the damping rate of the suspended particles (both quantities are defined more precisely below). 
In [3] we showed how the collision rate increases very rapidly when St exceeds a threshold value, 
due to fold caustics making the velocity field of the suspended particles multi- valued. In [3], 
which discussed initiation of rainfall from turbulent clouds, it was sufficient to use a single-scale 
flow model of the turbulent motion (described by a correlation length rj and correlation time r) 
because the Stokes number is never very large for particles suspended in terrestrial atmospheric 
clouds. However, in astrophysical contexts it is necessary to consider flows with large values of 
St, where the multi-scale aspect of turbulent flow [3] becomes important. (It is hard to study 
cases where St is large in terrestrial contexts because heavy particles fall out of the fluid) . 

In the following we derive an expression for the collision rate in a highly turbulent flow with 
large Stokes number. We employ the Kolmogorov cascade principle to deduce an expression for 
the variance of the relative velocities of colliding particles, which in turn determines the collision 
rate. 

Formulation of the problem. We assume that the drag force on a particle is proportional to 
the difference in velocity between the particle and the surrounding gas, so that the equation of 
motion is 



r = j[u(r, t) — r] 



(1) 



where r is the position of the particle and u(x,t) is the fluid velocity field (until the particles 
come into contact). This equation is familiar in the context of Stokes's law for the drag on a 
sphere, where the damping rate 7 is proportional to the kinematic viscosity v. In astrophysical 
applications, the mean free path of the gas is typically very large compared to the size of the 
particles [2], but equation (JTJ) remains applicable [5]. In this 'Epstein regime' the damping rate 
given by 7 = cp g / p p a, where a is the radius, c is the mean molecular speed of the gas and p g , 
p p are the densities of the gas and the particles respectively [5]. 

In [3] we demonstrated that the rate of collision R for a single suspended particle may be 
well approximated by 

R = RdiS + #adv + exp(-A/St)fl gas . (2) 

Here i?diff is a ra te of collision due to Brownian diffusion (and therefore independent of the 
intensity of the turbulence), i? a dv is the rate of collision due to the shearing effect of the flow, 
described by Saffman and Turner [6] and -R gas is the collision rate predicted by a 'gas-kinetic' 
model, introduced by Abrahamson [7], in which the suspended particles move with velocities 
which become uncorrelated with each other and with the gas flow. The exponential term de- 
scribes the fraction of the coordinate space for which the velocity field is multi-valued and A 
is a 'universal' dimensionless constant. The exponential dependence of the rate of caustic pro- 
duction on St was noted in [8] and recent simulations of Navier-Stokes turbulence suggest that 
A « 2 [9]. The rate R gas greatly exceeds -R a dv and Rdis, but the gas-kinetic theory is only 
applicable when the velocity field of the suspended particles is multi-valued. The mechanism for 
the particle velocity field becoming multi- valued is the formation of fold caustics, described in 
[3]. The formation of caustics can be modelled as a process of diffusion-driven escape from an 
basin of attraction, similar to the Kramers model for a chemical reaction [Sj. The exponential 
term is therefore analogous to the Arrhenius term exp(—E/kT) in the expression for the rate of 
an activated chemical reaction. The abrupt increase of the collision rate as the Stokes number 
exceeds a threshold value was first noted in numerical experiments by Sundaram and Collins 
[TP] . A theory proposed by Falkovich et al. [UJ emphasised the significance of caustics, but 
does not allow accurate quantitative comparisons with numerically calculated collision rates. 

Accounting for the multi-scale nature of the flow will change the expression for the gas- kinetic 
collision rate i? gas in ([2]) but will leave the others unchanged. For a gas with particle density no, 
this contribution to the rate of collision between particles of radius a is 

-Rgas = 4ira 2 n (Av) (3) 

where Av is the relative velocity of two suspended particles at the same position in space (and 
angular brackets denote averages throughout). For collisions in a conventional fluid a 'collision 
efficiency' factor is included to account for deflection of particles by the fluid trapped between 
them (see, e.g. [6]), but in the Epstein regime this factor is not required. Suspended particles 
exhibit a tendency to cluster when St ~ 1 [121 [T3] and in some circumstances the density no 
might have to be modified to take account of this effect 

Relative velocities in a multiscale flow. In the following rj and r are taken to be the dissipative 
length and timescales of the flow. According to the Kolmogorov theory of turbulence, these 
quantities are determined by the rate of dissipation per unit mass, £, and the kinematic viscosity 
v: we define 

' T= {s) ■ (4) 

If the turbulent motion is driven by forces acting on a lengthscale L ^> rj, the velocity fluctuations 
of the fluid have a power-law spectrum for wavenumbers between 1/L and I/77 [3]. 

In a multi-scale turbulent flow, when 7T <C 1 the motion of the suspended particles is 
underdamped relative to the motion on the dissipative scale, but (unless 7 is smaller than 
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frequency scale of the largest eddies) it is overdamped relative to slower long-wavelength motions 
in the fluid. The relative velocity of two nearby particles is a result of the different histories of 
the particles. If we follow the particles far back in time to when they had a large separation, their 
velocities were very different, but these velocity differences are damped out when the particles 
approach each other. A stochastic model of this process was discussed in |14| . which gave an 
expression for the relative velocity of colliding particles in terms of the velocity and length scales 
of the largest eddies. Here we show how to surmise the variance of the relative velocities by using 
the Kolmogorov cascade principle. Our result is universally applicable, in that it is expressed in 
terms of the rate of dissipation per unit mass, £, rather than in terms of the particular nature 
of the driving process. 

The motion of two suspended particles is determined by their damping rates 71, 72 and by 
the properties of the velocity field. When the particles are underdamped relative to the smallest 
dissipative scale, but overdamped relative to the 'integral' (driving) scale, we can apply the 
Kolmogorov principle [4] , that motion in the inertial range is independent of the mechanism of 
dissipation (i.e. it is determined by the rate of dissipation £ but it does not depend on v). The 
moments therefore depend only upon £ and 71 and 72. For the second moment of the relative 
velocity, dimensional consideration then imply 

(Av 2 ) oc £/j (5) 

where 7 must be some weighted average of 71 and 72, given by a formula which is symmetric 
under interchange of labels. 

We can gain some information about the form of the weighted average 7 by considering the 
relative velocities of particles with very different damping rates. In the case where one particle 
is much more heavily damped, 72/71 3> 1 say, the particles with damping rate 72 may be treated 
as if they are advected with the flow. Thus Av = \v\, where v = r — u is the velocity of the 
particle with damping rate 71 relative to the surrounding fluid. From equation ([1]), this satisfies 
v = u — 71 v, which has the solution 

v(t) = I dt' exp[7i(i'-t)]w(i') 

/t n„. i-t 

dt' exp[ 7 i(t' - t)] — (t') + / dt' exp[ 7l (t' - t)](v ■ V)u(r(t'),t') (6) 

where Du/Dt = du/dt + (u ■ V)tt is the Lagrangian acceleration of the fluid, which fluctuates 
on a timescale r. If \v\ ^> ij/t, the integrand of the final term in © fluctuates on a timescale 
r]/\v\ <C t, because in this limit fluctuations of u(r(t'),t') are determined by the rate of change 
of its first argument. Under this condition on the typical size of the velocity v (which is verified 
for 71T <C 1 below) the final integral in equation ([6]) may be neglected because its integrand 
fluctuates very rapidly. In the limit as 71 r — > 0, the variance of the velocity of the particle 
relative to the fluid approaches 

, 9 . 1 f°° I Du , . Du , . \ 

Noting that £t = (r//r) 2 , we confirm that (v 2 ) 3> (t]/t) 2 when 71 r <C 1. The dimensional 
arguments of Kolmogorov theory imply that I is a function of £ and v\ observing that X has 
the same dimensions as £ and we conclude that X oc £. Kolmogorov's 1941 theory of turbulence 
suggests that we should write X = K£ with K a universal constant, but in practice K may have 
a weak dependence upon Reynolds number due to intermittency effects [1]. Thus we have 

(A« 2 > = — = — (8) 
V 1 111 2 71 1 ; 
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when 72/71 3> 1 • Equations © and (jHJ) are consistent if we write 



( A «a> = -O ^(|ln( 7 i/72)|) (9) 

27172 

where g(x) is everywhere positive and approaches a finite limit as x — > 0. Furthermore 

lim 5 (x) = 1 (10) 

for consistency with (jHJ). 

In the case where the values of 7 are very different, we expect that the probability density 
function of the relative velocity is a three-dimensional Gaussian (or Maxwell-Boltzmann) dis- 
tribution, because the central-limit theorem is applicable to equation (jSJ) in the limit St — > 00. 
This implies that (Av) = ^8/37Ty/ (Av 2 ) . In the case where the damping rates are comparable, 
it would be of interest to evaluate the distribution function from a mechanistic model. 

Conclusions. Equations (JSJ) , (|2|) and fl3J) give the collision rate for underdamped particles in a 
highly turbulent flow, in terms of a measure of the turbulence intensity X = K£. The constant of 
proportionality K could be evaluated by direct numerical simulation of turbulent Navier-Stokes 
flow. The collision rate R is asymptotically proportional to \/E as the turbulence intensity £ is 
increased, provided that the Stokes number of the suspended particles is very large. The results 
will be of value in producing reliable estimates the rate of collision of dust particles in accretion 
discs. 
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